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Abstract 

We present an overview of the differential geometry of curves and surfaces using 

o. 

^ . examples from soft matter as illustrations. The presentation requires a background 

03 . 

only in vector calculus and is otherwise self-contained. 
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I. INTRODUCTION 

Though geometry is a common part of our early schooling, a rigorous and thorough 
education in physics usually tries to expunge it from our thought. Because of their predictive 
powers, there is good reason to emphasize numbers and formulae. Analytic geometry is the 
usual emphasis, while classical geometry is relegated to popular expositions. Differential 
geometry is a bridge between shapes and analytic expressions and is often the appropriate 
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language for modern physics. Nonetheless, when necessary, geometry is often slipped in as 
a bitter or, at least, tasteless pill - just enough is presented to get on with the analysis. In 
these lectures I have made an attempt to introduce the basics of differential geometry in the 
style of a mathematics text: the ideas are grouped by mathematical subject as opposed to 
physical subject. Nonetheless, I follow each newly developed topic with an example from 
soft matter, not only to illustrate the usefulness of the mathematical structure but also to 
aid the reader with a concrete example. 

This is by no means a textbook and many details of the physics are left for the diligent 
reader to find in the references. For technical details of many of the topics discussed here, 
the reader is referred to Elements of Differential Geometry by Millman and Parker flMillman 
and Parker, 1977|) , though any standard reference on classical differential geometry should 
suffice. A word on notation: I have tried throughout to explicitly display the functional 
dependence of all the fields and functions in formulae. However, sometimes this would make 
the notation awkward and the dependencies are dropped. In each case the context should 
make clear any lack of precision. 

Finally, I have tried to reduce as much as possible the use of the powerful formalism of 
differential geometry. While this formalism is useful for performing complex calculations 
unambiguously, great expertise is often required to extract the geometrical and physical 
meaning of these calculations. An excellent complement to these notes are the lecture 
notes by Frangois David ( David, 1989| ) which present the mathematical elegance and logical 
compactness of this subject. 



II. LOCAL THEORY OF CURVES 



A. Conformations of Polymers: Motivating a Geometrical Description of Curves 

Random walks abound in physics. They are the basis for understanding diffusion, heat 
flow and polymers. However, polymers are the most interesting of the three: polymers, unlike 
diffusing particles leave a "tail" behind them which they must avoid. They are described 
by self-avoiding random walks. As an introduction to the power of geometrical modelling, 
in this section we will consider the behavior of stiff polymers at the shortest lengthscale 
amenable to a continuum analysis. At these scales, polymers are not random walks at all, 
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but should be thought of as stiff rods. We take this as our starting point. When we describe 
a polymer as stiff, we are ascribing to it an energy cost for being bent. To model this, we 
consider a curve R(s), parameterized by its arclength s and construct the tangent vector 
to our curve T(s) = dR(s)/ds at a point s on the curve. In the next section we will show 
that the tangent vector is of unit length. If the tangent vector is constant along the curve 
then it is a straight line and does not bend. Thus the energy should depend on derivatives 
of the unit tangent vector T(s). Indeed we call the magnitude dT(s)/ds the curvature of 
the curve. In the next section we will discuss the geometry of curves in greater detail. In 
the meantime, we write the energy as: 



E - l A 

-^curv 2 



L 



dT(s) 



-i 2 



ds, (1) 



ds 

where s is a parameter that measures the arclength of the curve, and A is a measure of the 
stiffness. To study the statistical mechanics of the curve we write the partition function for 
T: 

f T(L)=T L 

Z (T L ) = / [ dX ] e -Bcu f v[T]/fc B T ( 2 ) 

JT(0)=2 

This is the partition function for the curve which starts with its tangent vector along the 
^-direction and ends with its tangent vector equal to T^. Note that if we can calculate 
( T(s)T(s') ), then we can integrate with respect to s and s' to obtain: 

Calculating the correlation function of the tangent vectors proves to be straightforward. 
Though there are many ways to proceed flDoi and Edwards, 1986|) , we will choose here an 
analogy with quantum mechanics. In ([|) let s = it. Then the partition function becomes: 

r T(L)=T L f • r -iL * / r 7 T \2 > ) 

Z(T " i) = / T( oM ldT] ^\^l dt l{n)\ (4) 

we recognize this as the path integral solution to Schrodinger's equation for a single particle 
where k B T replaces h, A replaces the particle's mass, and T is the position of the particle. 
Since T(t) lives on the unit sphere, this is just quantum mechanics on a sphere. The 
Schrodinger equation is: 

KTdz (k B T) 2 f2r7 

~—dL = —2A~ L Z > (5) 
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where L is the angular momentum operator. Changing back to our original coordinate s 
and defining L p = A/k B T we have 

dZ{T,s) 1 



ds 2L P 



L*Z(T,s). (6) 



We will soon discover that L p has an important interpretation. 

Since the polymer is the same all along its length, we have ( T(s) • T(s') ) = ( T(s — 
s') ■ T(0) ). Using polar coordinates and enforcing the limits of integration in @ so that 
T(0) = z, we discover that we are interested in (cos 8(S) ) at S = s — s': 

f_, g?(cos 0) cos #Z(cos S) 

(cos9(S)} = 1 (7) 

J i rf(cos 6)Z(cos 6, S) 

where we have been sure to normalize Z to get a probability. The expectation value we 
seek satisfies a differential equation since Z(cos9, S) satisfies (^). Moreover, since L 2 is a 
Hermitian operator, we have 

d -^mi = J r( e cosl){ S)) = -U C os l >(S)), (8) 
where we have used the fact that L 2 cos^ = —2 cos# and L 2 (l) = 0. It follows that 

(cos[0(s)-0(s')]> =e- | - i 'l /I ». (9) 
Integrating this expression as in (Q), we have 

(R(L) - R(0)) 2 ) = 2L P (L-L p + L p e- L/L ?) . (10) 



From ([jp and ( PD[ ) we glean the meaning of L p . The first equation shows that the tangent 
vectors along the curve are uncorrelated after a distance L p . For this reason L p is called 
the persistence length ( |de Gennes, 1970] ). The second equation shows us that for L much 
longer than L p , a stiff rod behaves as a random walk: i.e. the average square distance that 
is travelled is proportional to the number of steps or length of the walk, R 2 oc LL P . For L 
much shorter than L p , we may expand (|ToD to see that R 2 oc L 2 . 

The physics of stiff rods can be used to study other phenomena. For instance, vortices in 
fluids, superfluids and superconductors have a bending stiffness arising from self-interactions. 
In the past decade the physics of stiff rods has been adapted to study the mechanical 
properties of DNA ( Marko and Siggia, 1994| , |1995| ), and has been augmented to include 



FIG. 1 A curve in space. At R there is a tangent vector T, a normal vector N and the binormal 
vector B. 



twisting degrees of freedom ([Kamien et al, 1997j [Marko, 1997| ), the topological constraints 
imposed by self- avoidance flVologodskii et al, 1992| ) and by closed loops flMoroz and Kamien 



|1997| ; |Moroz and Nelson, 1997] , |1998|) , and the effects of sequence disorder ( Pensimon et al. 
19981; INelson, 1998;). 



B. The Frenet-Serret Apparatus: DNA and Other Chiral Polymers 



In the previous section we have seen that a simple geometrical description of polymers 
leads to a precise description of their conformational behavior in a variety of regimes. It is 
now time to discuss more carefully the geometry that went into making the expression for en- 
ergy in (fj). A curve in three dimensions is a vector- valued function R(s) = [x(s), y(s), z(s)] 
that depends on s, a parameter that runs along the curve. Though we may choose to label 
points along the curve as we wish, it is usually most convenient to let s be the arclength 
along the curve. We will see how this simplifies our equations shortly. If the curve is L long, 
then s runs from to L. The first thing to do is to construct the unit tangent vector to the 
curve: as we have already asserted, the magnitude of the rate of change of the unit tangent 
vector is the curvature k of the curve. We have 

T(s) = iBr 
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where X'(s) denotes the derivative of X(s) with respect to s, and ||X|| is the length of the 
vector X. We have divided by the length of R'(s) to make T(s) a unit vector. However, 
if s measures the arclength along the curve then ||R'(s)|| = 1! To see this, recall that the 
length of a curve R(i) from t = to t — tf can be found by adding (integrating) the length 
of pieces of the curve together, each y/<iR(i) ■ d~R(t) long: 



If we choose t = s to be the arclength then we have 

-L 



ds\\R'(s)\\, (13) 

o 

where the upper limit is the same as the length of the curve. Differentiating both sides of 
(0), we see that ||R'(s)|| = 1. 

Having constructed the unit tangent vector, we may now take its derivatives to obtain 
the curvature. Since the derivative of a vector is another vector, we write 

T'(s) = «(s)N(s) (14) 

where k(s) is the curvature at s and N(s) is a new vector, the unit normal vector, which 
is also unit length, so that ||T'(s)|| = |ac(s)|. It is convention to choose k(s) to be always 
positive - the sign can alway be absorbed in the direction of N(s). Note that T(s) -N(s) = 
since the derivative of any unit vector is perpendicular to itself; if ||X(s)|| = 1 then 

lpC(.).X(.)] = l[l] 

2X(s) -X'(s) = 0. (15) 

Of course, as we move along s, the normal vector changes direction as well. Changes in the 
direction of the normal vector can come from two contributions: the normal can change by 
rotating towards or away from the tangent vector (of course the pair rotate together to stay 
perpendicular). It can also change by rotating around the tangent vector. The former case 
corresponds to the curve staying in the same, flat plane, while the second corresponds to 
rotations of the plane in which the curve lies at s. This plane is known as the osculating 
plane, from the Greek word for kissing. Moreover, since N'(s) is perpendicular to N(s), we 
must introduce a third unit vector to account for changes in the osculating plane. We choose 
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B(s) = T(s) x N(s) where x denotes the cross product, as shown in Figure 1. This vector, 
the binomial vector, is a unit vector perpendicular to both T(s) and N(s). We then have 



N'(s) = a(s)T(s) +t(s)B(s) 



(16) 



where a(s) is some function of s and r(s) is called the torsion of the curve. It is a measure 
of the rate of change of the osculating plane. Why don't we give a(s) a name? Because we 
note that by differentiating the relation T(s) • N(s) = we get: 

= T'(s) • N(s) + T(s) • N'(s) = k(s) + a(s), (17) 

so a(s) = —k(s) and ( p~6|) becomes: 

N'(s) = -«(s)T(s)+r(s)B(s). (18) 

Finally, we may calculate B'(s) to complete our analysis of the curve. We have 

B'(s) = T'(s) x N(s) + T(s) x N'(s) 

= k(s)N(s) x N(s) - k(s)T(s) x T(s) + r(s)T(s) x B(s) 

= -r(s)N(s) (19) 

where the last line follows from the rule ax (b x c) = b (a • c) — c (a • b) . Putting ( [Tip . 
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T( S )" 






N(a) 


(20) 




B( S ) 





p!8|), and (|T^) together we have the Frenet-Serret equations for a curve in three-dimensions: 



d_ 
ds 



This shows that k(s) is the rate of rotation of T(s) about B(s) and similarly, r(s) is the 
rate of rotation of N(s) about T(s). Written as one matrix equation, the Frenet-Serret 
formula tell us something very important: given a curvature k(s) and a torsion r(s), we 
can reconstruct our entire curve up to a translation (since we can change the origin) and 
a rotation (since we can rotate the initial orthonormal triad {T(0), N(0), B(0)}. Once we 
have set the location and orientation of that triad, however, the entire curve is determined 
by only two parameters, not three as one might have thought. 

There is a difficulty with the Frenet-Serret frame: when the curvature vanishes, N(s) 
is not well-defined and therefore B(s) and, more importantly, the torsion are not defined 
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either. Thus curves that have straight segments are problematic from the point of view of 
the Frenet-Serret frame. Moreover, if we consider a helix 

(- 



qs \ QS 
r cos I — , , r sin 



V^r) 2 + TJ ' VVW + V #F+1 
then the curvature and torsion are constant: 



(21) 



g 2 r 

Ms 



1 + (gr) 2 

• (22) 

As r — > we approach a straight line. However, though the curvature vanishes in this limit, 
the torsion does not! This is a problematic feature of the Frenet-Serret frame - since the 
torsion is the magnitude of the derivative of N(s), it is only a meaningful quantity when 
N(s) = T'(s)/k(s) is unambiguously defined, and this requires that k(s) ^ 0. In §IIIB, we 
will offer a different frame that does not suffer from this problem. 

Knowing that the there are only two parameters needed to describe a space curve, we 
can now augment (0) to include other effects. One interesting effect ( Harris et ai, 1999 ) is 



the behavior of chiral polymers. While the curvature does not distinguish between left and 



right handed helices, we can see from (22) that the torsion is sensitive to the sign of q. Thus 
we might add terms to ([!]) to favor a particular chirality for the stiff polymer. We might be 
tempted to write 

E = E curv + E* = jf ds j^ 2 (s) + | [r{s) - T ] 2 j . (23) 

Though this energy appears acceptable, and favors an average torsion r = r , it only ac- 
centuates the ambiguity of r when k = 0. Moreover, since (|IJ) is a functional of R(s), we 
should construct the new energy only in terms of R(s) and its derivatives. In addition to 
[R"(s)] 2 = [T'(s)] 2 , we can also construct the term (Kamien et gl, 19971 ; IMarko, 19971) 



E* = -* 
2 



/ dsT(s).[V(s)xT»(s)]. (24) 



From the Frenet-Serret formula, we then have: 



E = £ds |^ 2 (s) + ^k\s) - |k 2 ( S )t( S ) + ^ 2 (,)r 2 ( S )| , (25) 

where A, C, a and (3 are all positive. Note that in the special case of constant curvature 
and torsion, [T(s) x T"(s)] 2 = k 2 t 2 . We have added this term and the extra quartic term 



FIG. 2 The tangent spherical map. The tangent vector of the curve T traces out a curve T on the 
unit sphere. The length of T is the integrated curvature along the original curve. 

to E for reasons that will become clear in the following. This form does not suffer from 
the torsion ambiguity: when k(s) = the torsion drops out of (|25|). One can minimize 
this energy for the helix (pip , which has constant curvature and torsion, to find the ground 
state conformation of a chiral polymer. For appropriate parameter values, the tendency for 
torsion can overcome the tendency to be straight and both k and r will be nonvanishing. 

III. GLOBAL THEORY OF CURVES 

A. Fenchel's Theorem: Energetic Bounds on Closed Curves and Knots 

When we consider the energetics of closed curves, it is clear that the curvature cannot 
vanish everywhere (or the curve won't close on itself). We can actually establish a lower 
bound on the curvature energy by appealing to a theorem on the total curvature of the 
curve. Not only will this be useful when considering closed polymer loops, the technology 
of the proof will help us later on when we consider the geometry of surfaces. 

Fenchel's theorem states that the total curvature of a closed curve is at least 2tt. This 
is not unreasonable: a circle of radius R has curvature 1/R and so the total curvature is 
2itR/R = 2tt. In general, other closed curves only have more curvature (in fact, only planar 
convex curves have an integrated curvature of exactly 2n). The statement of the theorem 
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is: 



j> K(s)ds > 2tt (26) 



for any closed curve, where, as usual, s is the arclength of the curve. Before we prove this 
fact, note that via the Cauchy-Schwarz inequality we have 

jf (l)K(s)ds < Sj^ K 2 (s)dssjj^ l 2 ds. (27) 

Squaring both sides and using (^6|) we find that the curvature energy (0) satisfies 

1 C L L 

E CUTV = -k B TL p <h K 2 {s)ds > 2n 2 k B T-^. (28) 

2 Jo L 



To prove Fenchel's theorem we introduce the tangent spherical image ( [Millman and 



Parker, 1977[) of the curve: we take the tail of the unit tangent vector T(s) at R(s) and map 



it to the center of the unit sphere. As shown in Figure 2, the tip of the tangent vector then 
traces out a curve T on the surface of the unit sphere as s goes from to L. The differential 
of length of the curve on the tangent sphere is ||T(s)'||ds, in analogy with (|13|) . But the 
magnitude of T' is just the curvature, so we find that the length of T on the unit sphere is 

-L 

K(s)ds, (29) 



the total curvature of our original curve R(s). Now we note that by definition of T(s), 

R(s) - R(0) = / T(s')ds', (30) 
Jo 

so if the curve is closed, R(L) = R(0), and the left-hand side of ( |30"D is 0. This is really 
three equations, one for each component of T. It tells us that the curve T on the tangent 
sphere can never be in only one hemisphere. In order for the curve to close, it must turn 
around. But if the tangent spherical map were in one hemisphere, we could take it to be the 
upper hemisphere where T z (s) > 0. Yet if T z (s) is always positive then its integral cannot 
vanish! Therefore, given any hemisphere on the tangent sphere, the curve must be in it and 
its complement. But this means that T must be at least the length of the equator, 2ir. We 



have thus proven (26). The value of this result is not just the bound on the energetics of a 
closed polymer. It has introduced us to the tangent map. As we will see in the following it 
is very useful to take vectors off of curves and surfaces and translate them to the center of 
the unit sphere. This was our first taste of this procedure. 
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In closing this section we mention the Fary-Milnor theorem which pertains to the inte- 
grated curvature of a non-self- intersecting closed knot. Not surprisingly, this theorem states 
that a knot has to go around at least twice: 

-L 



o 



n{s)ds > 4tt. (31) 
Using the same reasoning that led to (^8|), we have 

r L L 

k B TL p f K 2 (s)ds > (47i) 2 k B T-f- (32) 
Jo L 

for a knotted closed curve. 

B. The Mermin-Ho Relation: Basis Vectors to the Rescue 

Up until this point, we have been able to study the geometry of lines without the in- 
troduction of a frame: that is, all our results rely on the original tangent vector T(s) and 
its derivatives - at no point were we required to choose a basis or coordinate system to 
calculate any quantities. Unfortunately, we cannot continue our discussion of lines and we 
certainly cannot discuss surfaces without introducing some more paraphernalia. We need 
to introduce a set of spatially varying basis vectors in order to define quantities in addition 
to the tangent vector, the curvature and the torsion. We need to make sure, however, that 
all of our physical quantities do not depend on our arbitrary choice. In this section we will 
derive and explain the Mermin-Ho relation (|Mermin and Ho, 1976|), originally derived in the 



context of superfluid 3 He-A, a phase characterized by an order parameter with two direc- 
tions, a "long" direction and a "short" direction perpendicular to it. Though this section 
is the most technical, it is also straightforward. Fortunately, once we have established this 
result we will be able to make use of it again and again in the following. 

Often we have a vector n(x) defined everywhere on a surface or in space. It could be 
the director of a liquid crystal, the normal to a surface or some other field of interest. We 
now are interested in a vector N(x) which is always perpendicular to n(x). This vector 
might point to the nearest neighbor or in some special direction on the surface. To define 
this vector, we introduce two new unit vectors, ei(x) and e 2 (x), which are defined to be 
perpendicular to n(x) and each other so that {ei(x), e 2 (x), n(x)} is an orthonormal triad 
and ei(x) x e 2 (x) = n(x). A vector which is defined relative to the spatially varying plane 
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defined by ei(x) and e2(x) will change not only because its true direction changes (relative 
to a fixed triad) but also because the basis vectors change. If N(x) is a unit vector then 
N(x) = cos[6'(x)]e 1 (x) + sin[0(x)]e 2 (x) is always perpendicular to n. Its derivatives in the % 
direction are vectors as well. In the plane normal to n their components are: 

ei(x) • ^N(x) = -sin0(x) [^0(x) - ei (x) ■ ^e 2 (x)] (33) 
e 2 (x) ■ <9iN(x) = cos0(x) [d;0(x) + e 2 (x) ■ o^x)] 

= cos0(x) [$0(x) - ei (x) • <9,e 2 (x)] (34) 

where we have used the fact that ei(x) • e 2 (x) = in the final equality. Notice that the 
derivatives of N(x) depend on both gradients of 0(x) as well as derivatives of the spatially 
varying basis vectors e$(x). Since the basis vectors were chosen arbitrarily, one might be 
concerned that the derivatives of N(x) are poorly defined. However, there is a concomitant 
change in 0(x) whenever the basis vectors change so that the gradients of N(x) are well 
defined. The problem is with gradients of 0(x). 

To disentangle the arbitrary dependence on basis vectors, we start by considering a vector 
N (x) which is constant in the instantaneous plane perpendicular to n(x) so that (|33|) and 
([34]) both vanish. For N D (x) to be constant in the plane perpendicular to n(x), 0(x) must 
be equal to some o (x) with V6 l Q (x) = ei(x) • Ve 2 (x) = f2(x), where the last equality 
defines a new vector field fi(x) = e x (x) ■ Ve 2 (x) called the spin connection. Note that 
we can only solve this equation for o (x) if V x Q = 0. If we now consider an arbitary 
vector perpendicular to n(x), N(x), we should focus not on gradients of 0(x), but rather 
on [0(x) — 8 (x)}, which is a measure of how much 0(x) deviates from its "constant" value. 
But V0 o (x) = f2(x), so V [0(x) - o (x)] = V0(x) - f2(x). Thus we see that by subtracting 
f2(x) from V0(x), we remove that part of 0(x) that is induced by a spatially varying basis. 
Moreover, even if V x f2 is nonvanishing, we can generalize this discussion to form the 
combination D0(x) = V0(x) — fi(x), the covariant derivative. It measures the true changes 
in 0(x), relative to o (x). If [0(x) - o (x)] is a smooth field then V x V [0(x) - o (x)] = 0, 
or 

V x V0(x) = V x n(x). (35) 

We might be tempted to set the left hand side equal to zero, since usually the curl of a 
gradient vanishes. However, as we are about to show, the curl of fi(x) is not always zero! 
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This means that there must be some sort of singularities or defects in #(x) (and # (x) - the 
difference #(x) — 8 (x) is smooth). We will talk more about defects in §Vb. 
We have for the i th component of the curl: 

[V x n(x)] i = e l]k d 3 [e?(x)0 fc e£(x)] 

= e ijk [dje™ (x)] [<9 fc e£ (x)] + e? (x-)e ijk djd k e% (x) 

= e ufc [^e?(x)][^(x)] + (36) 

where the last term vanishes due to the antisymmetry of and we have used indices to 
make the calculation unambiguous 1 . We have used the Einstein summation convention of 
dropping summation signs for repeated indices. Unless otherwise indicated, an index should 
be summed over if it appears twice in any formula. Both the Greek and Roman indices run 
from 1 to 3. Now consider the object o^e^x). Since ei(x) is a unit vector, its derivative is 
perpendicular to it. Therefore, we can write c^e" (x) in terms of the basis vectors n a (x) and 
4(x): 

a,-e?(x) = A,(xK(x) + fl,-(x)e£(x). (37) 
Similarly, we can do the same for dj.e% (x): 

fc e£(x) = C fc (x)n a (x) + D fc (x)e?(x). (38) 

Putting these expressions into (|36|) and using the orthogonality of our triad, we find: 

[VxO(x)] ?; = e l , fe A J (x)C fc (x). (39) 

By taking the dot product of equation (j37|) with n(x), we find A,(x) = n /3 (x)9 ? e^(x) = 
— (x)(9 J n /3 (x). Similarly, C/-(x) = — e2(x)<9fcn 7 (x) and so 

[V x fl(x)] < = 4 (x)^(x)cy fc [a,-n"(x)] [^(x)] . (40) 



1 The only tensor on which we will rely is the antisymmetric tensor eij k - It is defined by 

T if ijk is an even permutation of 123 
Cijk — ^ — 1 if ijk is an odd permutation of 123 

if any two of i,j or k are the the same 
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We see then that the curl of fi(x) does not always vanish, but it appears that it depends 
on our arbitrary vectors ei(x) and e2(x). Note however that if we interchange the indices (3 
and 7 in (f40|) that it is equivalent to interchanging j and k, which would introduce a minus 
sign. Thus we have 

[V x 0(x)], t = ~ [e?(x)eZ(x) - eftxjef (x)] e ijk [«V(x)] [^(x)] 

= ^e a/3l n a (x.)e ijk d j n l3 (x)d k n' y (x), (41) 

where we have used the orthonormality of {e^x), e 2 (x), n(x)}. This is the celebrated 
Mermin-Ho relation flMermin and Ho, 1976|) . Note that because of the antisymmetry of 



€ij k , it is unnecessary to keep the brackets around the gradients of n(x). 

This relation between n(x) and vectors perpendicular to n(x) is rather remarkable. Note 
that our choice of e^x) and e 2 (x) is arbitrary so that fi(x) is not a constant of the geometry. 
However, ([41]) shows that VxO(x) only depends on n(x), and not our choice of basis vectors! 
This is why the Mermin-Ho relation has been introduced in this section on global properties 
( [Kamien, 1998| ): if we have a closed curve T of length L, then 



/ [V0(x) - fi(x)] • dR = [6(L) - 9(0)] - [[ [V x n(x)] • dS 
Jr JJm 

= [6(L) - 9(0)] - J J U, up e aPl n a (^)d u n^(^)dy(^)dS^2) 

where we have used Stokes theorem to change an integral around a curve V into an integral 
over a capping surface M. In fl4*2|) dS is an element of area. Thus if we consider changes in 
a vector around a closed curve, our choice of ej(x) is irrelevant. Though it is necessary that 
N(L) = N(0), 9(L) — 9(0) can change by an integral multiple of 2ir. The rest of the change 
in the direction of N(x) comes from the geometry of n(x). We will see how this extra term 
plays an important role in the next section, in the physics of surfaces and the physics of 
defects in three dimensions. 

C. Link, Twist and Writhe: Dynamics of twist-storing polymers 

One of the more interesting stiff polymers is DNA. Though it is well known that it has 
great biological significance ( |Watson and Crick, 1953 ) , it is of interest in materials for at 



least two other reasons. First, the persistence length of DNA is quite long, roughly 50 nm. 
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More importantly, because there are a plethora of enzymes available to cleave DNA, a sample 
of monodisperse polymers can be prepared much more readily than in any synthetic system. 
We have already discussed conformations of chiral polymers like DNA in §11. However DNA 
has one more interesting element: it can form into a closed loop. This is interesting because 
DNA is actually a double stranded helix. Therefore the number of times that one strand 
wraps around the other is fixed when the loop is closed (when the strand is open the helix 
can unravel). Thus there is a conserved quantity and this leads to a constraint on the 
possible dynamics of the double strand. More generally, any polymers that cannot unwind 
along their axis are known as twist-storing polymers. 

When we have two closed curves T and T', we can assign a linking number Lk to them 
which counts the number of times one loop passes through the other. There is a simple way 
to calculate this given the two curves R(s) and R'(s) using an analogy with Ampere's law. 
If we think of the first curve T as being a wire carrying a current /, then we know that the 
line integral of the generated magnetic field B around the closed curve T' is Arcnl, where n 
is the number of times that the current passes through the closed loop 2 . Setting / = 1, we 
have 

Lk = n=^- I B(x') • dx. (43) 



4ir JT , 

We may calculate the resulting field B(x') from the wire by use of the Biot-Savart Law: 

A fldixr jfdxxfx'-x] .... 

B ^ = f r -^ = f r |„_ X |. • < 44 > 

Putting these together we find that (with the curves having length L and L', respectively): 

1 / / [x-x']' (rfx x dx 1 ) 



4"7T J r J T i |x — x'| 3 

fL' 

ds' 



dR(s) dR'(s' 

X 



[Rjf) - RV)j 

>(a'\\3 



R(s) -R'(s 



= G(T,T') (45) 

the last equality defines G(T,V), the Gauss Invariant. Thus two closed curves that cannot 
separate and rejoin must keep this invariant constant. 

While this may be elegant, it is not especially useful when studying molecules like DNA. 
At long lengthscales, DNA appears as a single filament. It would be useful to recast the 



For simplicity we work in cgs units with c = 1 
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linking number in terms of the single polymer picture of DNA. To do this, we consider two 
curves that are close together ([Fuller, 1971]) . The first curve is R(s), while the second curve 
is 



R'(s') = R(s') + eu(s') 



(46) 



where, as is the tradition, e is a small number and u(s') is a unit vector that is perpendicular 
to the tangent vector dH(s')/ds = T(s'). Inserting these expressions into (EH), we have 



Lk= — I ds I ds'T(s) 
4vr J J 



T(s') + 6 



etuis' 
ds 



[R(s) - R(s') - eu(s')] 
iRfs) -R(s') -eu(s')l 3 ' 



(47) 



We now want to take e to zero. While the numerator in (H7T) does not make this a problem, 
note that the denominator diverges whenever s = s'. Thus as long as \s — s'\ > 5 we can 
take e -> flKleinert, 199CQ . We thus have 



Li 



47T 



(is' 



" S ds+ r J T W xTM-[RW-RM] 

+5 



iRfs) -R(s')l 3 



+— ds 7 

47T 



c/sTfs) x 



T(a') + e 



du(s') 
ds 



\R(s) -R(s') 



euis 



|R(s) - R(s') - eu(s') 



(48) 



The first integral bears a resemblance to the Gauss invariant, while the second integral 
depends on the vector u(s). Since S is small, we can expand R(s) around s' to find: 



and 



R(s) = R(s') + (s- s')T(s') + 



dT(s') 

T(s) = T(s') + (s - s')^fl + 

ds 



(49) 



(50) 



Inserting this into the second integrand and using the fact that T(s) • u(s) = 0, we find 

Lk = — t ds t ds'T(s) x T(s') ■ R (*) ~ R ( g/ ) 
4tt7 X 1 ; 1 J |R( S )-R( S ')| 3 



ds' / ds 

47T j 



s'+S 



T(s') +(s-s> 



dT(s') 
ds 



T(s')+e 



du(s') 
ds 



ems 



(51) 



[{s - s') 2 + e 2 ] 3/2 

where we have taken e — > in the first integral and have used T(s') x T(s') = in the 
second. Note that the s dependence of the second integral is manifest and we may do the 
integration over s, let e go to zero first and then let 5 go to zero. We have separated the 
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FIG. 3 The tangent spherical map. The area of the patch M on the sphere is, through Stokes 
theorem, the line integral around T of ei • Ve2- 



link into two integrals. The first is called the writhe, Wr, and it only depends on the curve 
R(s). The second is called the twist, Tw, and it depends on u(s). We have 



Wr 



Tw 



dsi ds'T(s) xT(s') ■ 



Air 



— <b dsT(s) ■ 



R(s) - R(s') 
K(s) -R(s'W 



du(s) 



US X 



ds 



(52) 
(53) 



and thus we arrive at Fuller's celebrated result ( |Fuller, 1971| ): 



Lk=Tw+Wr. (54) 

Though the expression for writhe ( |52"D bares a strong resemblance to that for the Gauss 
invariant ([15]), it is not the same. In the Gauss invariant we were considering two different 
curves that did not touch and there was no need to expunge a singularity from the integra- 
tion. The expression for writhe, on the other hand, is for the same curve and is only defined 
in the limit described above. Note further that the writhe is nonlocal, while the twist is 
local. This is the price we must pay: the link is topological, the writhe depends only on the 
backbone R(s) but is nonlocal, and the twist is local but we must know about u(s), i.e. the 
other curve. 



Note that the expression for the twist (^) suggests a use of the results of §IIIB. Since 
u(s) is perpendicular to T(s), we may expand it in a set of basis vectors e 1 (s) and e 2 (s) 
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both perpendicular to T(s), so that u(s) = cos#(s)ei(s) + sin 6*(s)e 2 (s). Then 



T s 



USX 



d s 0(s) - ei(s) • d s e 2 (s) (55) 



where we have used the fact that ei x e 2 = T (plus cyclic permutations). But this is the 
form we considered in the discussion of the change in a vector. Following the argument at 
the end of the previous section, we find that 

Tw = / ds [d s 9(s) - ei(s) • <9 s e 2 (s)] 

= m- J J ^ e ^ p e^ 7 T Q (x)^T^(x)^(x)^ 

= Lk- Wr, (56) 

where m is an integer (since 9{L) — 6(0) = 2irm). We would be tempted to identify the 
second integral Wr as the writhe, since we have an equation that reads justs as ([54]). This 
is not always correct and depends on our choice of basis vectors. However, it does establish 
a weaker result of Fuller's ( [Fuller, 1978|) that the writhe is the last term in (|56"D mod 1. As 



we will see in the next section, the integral has a geometric meaning: it is the area swept 
out by the tangent curve on the tangent spherical map, as shown in Figure 3. 

The difficulty with writhe is that it is non-local and therefore cannot be easily included 
in a local set of dynamical laws - the conformation of the entire curve must be known to 



calculate writhe. Moreover, writhe is only defined for closed curves and so the identity (p4|) 
would not appear to apply to open strands. While the total amount of writhe including 
the integral part is important for calculating ground states of a particular twisted ribbon 



( pain and Rudnick, 1999[ [Fain et ai, 1996] ; [Julicher, 1994j ; Marko and Siggia, 1994|) when 



considering changes in writhe, the constant integer part is less important. If the timescale 
for the diffusion of twist along the polymer is long enough, then even an open strand should 
feel the constraint of conserved link. 

To this end, a useful result for the change in writhe of a curve as a function of time t 
( lAldinger et aL, 1995|) : 

d t Wr(t) = i- IdsT(s,t)-[d t T(s,t) x d s T(s,t)] , (57) 

can be used. Note that this result follows from the preceding discussion: if we choose a 
coordinate system, then ( |5"7P (multiplied by dt) is the differential of area swept out by the 
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tangents of two closed curves with tangent vectors T(s, 0) and T(s, dt). Since the twist is 
really the local torsional strain of the polymer, we denote the twist density as u(s,t) and 
then we have 

d t Lk= jds {d t u(s,t) +d t T(s,t)-[d s T(s,t) x T(s,t))} + d t n (58) 

where n is the integer difference between Wr and Wr. Since continuous evolution cannot 
lead to discontinuous changes in the integer n and since link is conserved we have 

= jds {d t u + d t T-[d 3 T x T]}. (59) 

This conservation law need not be satisfied locally: the curve can twist in one place and 
writhe at some distant location to satisfy fl5P|) . In addition, the integer can change if the 
curve develops cusps and evolves in a non-smooth way. Because physics is local, however, 
we might expect that linking number is locally conserved and changes via a "link current" 
j. This would lead to a local conservation law 

d a j = d t u + d t T-[d B TxT\, (60) 

which satisfies (j59|) . This conservation law enforces total link conservation since § ds d s j = 
0, yet allows for local deviations in the twist and writhe ([Kamien, 1998| ) . 

This geometrically inspired conservation law has been verified in numerical experiments 
( |Goldstein et gi, 1998| |Wolgemuth et al, 2000|) which show that there are two modes of 
relaxation in the dissipative dynamics: twirling, in which the link relaxes through torsional 
modes along the polymer (like a speedometer cable), and whirling, in which the link relaxes 
through crank-shaft like motions of the entire chain. 

IV. LOCAL THEORY OF SURFACES 
A. The Area Element: Minimal Surfaces 

While it is easy to generalize the equation of a curve R(s) to an equation for a surface 
X(si,S2), it turns out that the theory of surfaces is much richer than that for curves. In 
the first place, there is no way to talk about arclength when defining s\ and S2, so we will 



20 




FIG. 4 The magnitude of the cross product of two vectors is the area of the parallelogram swept 
out by them. 

take them to be arbitrary coordinates of the surface. We can, however, consider the area 
of a small patch of our surface. To do this, we need to construct vectors tangent to the 
surface. As shown in Figure 4, if we have two vectors a(si, s 2 ) and b(sx, s 2 ) 111 the surface, 
then they span a surface area AS = |a(si, s 2 )| |b(sj, s 2 )| sin#, where 9 is the angle between 
the two vectors. We recognize this expression for AS* as the magnitude of the cross product 
a(si, s 2 ) xb(si, s 2 ). Since both of these vectors are tangent to the surface, their cross product 
is parallel to the unit layer normal n(sj, S2). Thus, we have 



AS = |n(>i, s 2 ) • [a(si, s 2 ) x b(si, s 2 ) 



(61) 



But, with this construction we can also construct the unit layer normal out of a(si, s 2 ) and 
b(si,s 2 ): 

a(si,s 2 ) x b(si,s 2 ) 



n(Si,s 2 



|a(si, s 2 ) x b(si, s 2 ) 



Putting this expression together with (|H]), we have 

[a(si, s 2 ) x b(si, s 2 )] ■ [a(si, s 2 ) x b(s a , s 2 )] 



|a(si, s 2 ) x b(si, s 2 ) 



[a(si, s 2 ) x b(si, s 2 



^ a 2 (si, s 2 )b 2 (si, s 2 ) - (a(si, s 2 ) • b(si, s 2 ))' 



(62) 



(63) 



We now need to construct two, nonparallel vectors a(si,s 2 ) and b(si,s 2 ). The surface 
provides us with two: a(si, s 2 ) = 9 Sl X(si, s 2 )dsi and b(si, s 2 ) = 9 S2 X(si, s 2 )c?S2. This tells 
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us that the differential area element is; 



dS = y (<9 sl X(si,s 2 )) (9 aa X(si, s 2 )) - (9 ai X(si,s 2 ) •9 a3 X(si,s 2 )) dsids 2 . (64) 

Often we also need the vector surface element c?S = n(sj, s 2 )dS, for instance when calculating 
electric and magnetic flux. The area of the whole surface M is simply 



A= / J {d Sl X) 2 (d S2 X) 2 - («9 Sl X • d S2 X) 2 d Sl ds, 



(65) 



What happens if we choose to parameterize the surface in terms of a new set of coordinates? 
If we have two new coordinates o\ and cr 2 defined by crj = Oi(si, s 2 ), then by the chain rule 

<9X <9X 9cr,- 9cr,- 



<9,,X 



J_ _ Z_J_Q X 

(9sj (9(Tj 9sj dsi ° 3 



(66) 



where, as usual, there is an implicit sum over j. Because of this sum, it appears that it 
would be somewhat tedious to reexpress ( |65| ) in terms of the new coordinates. However, we 
note that if we define a matrix g: 



<9 Sl X-<9 Sl X <9 Sl X-<9 S2 X 



<9 S2 x-<9 sl x <9 S2 x-d S2 x 



(67) 



then the expression in the radical of ( j65| ) is simply the determinant of this matrix. In 
component form, this matrix is gij(si, s 2 ) = d Si X.(si, s 2 ) • £? a .X(si, s 2 ). This matrix is known 
as the metric tensor which some people find useful ( |Misner et al., 19"73|) . Moreover, we note 
that the transformation ( |B"6] ) amounts to matrix multiplication of g(si, s 2 ): 



Qkmipx, ^2), 



(68) 



where gkm(&i, 02) is the corresponding matrix in the new coordinate system. Defining a new 
matrix O by Oij = then g = T gO. Thus when we change coordinates, we have 

A = / a/ det g ds 1 ds 2 = / a/ det T gO dsids 2 = / A/detg |det 0| dsids 2 . (69) 
Jm ./m 7m 

But the determinant of O is just the Jacobian determinant of the transformation from s< to 

a,, and so |det 0| dsids 2 = do\da 2 or 



A= a/ det g ds\ds 2 = / a/ det g da\da 2 

JM JM 



(70) 
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and so the area is invariant under coordinate transformations. This invariance (known 
as diffeomorphism covariance by the cognoscenti) is useful as it allows us to choose the 
most convenient set of coordinates for a given problem. One often used choice is the so- 
called Monge gauge or height representation, where S\ and S2 are chosen to be the x and y 
components of the surface vector X(s 1; s 2 ) an d the ^-component is a function of x and y: 

x = Xi(sx,s 2 ) = Si 
y = X 2 {sx,s 2 ) = s 2 

z = X 3 (sx,s 2 ) = h(sx,s 2 ) = h(x,y) (71) 

It is straightforward to check that the expression for the area (|65D becomes the more familiar 
expression: 



dh\ 2 fdh^ 2 



A -!J 1+ \ji) + W My (72) 

The height representation requires that for every x and y, there is only one height or, in 
other words, the surface can have no overlaps. Sometimes it is necessary to break the surface 
up into different regions or patches in order to make a good height representation. 

The height representation proves useful when considering the energetics of fluid mem- 
branes. By fluid we mean that the membrane has no internal structure and the molecules 
in it can flow freely. Often we are interested in the shape and dynamics of these membranes 
when they are under a uniform applied tension. Soap films are a common example, though 
in many cases the interface between two fluids or two distinct phases can also be thought 
of as a membrane under tension. The free energy of this system is simply F = $A, where 
<j is the surface tension and A is the area of the interface. Minimizing the free energy thus 
amounts to minimizing the area of the surface. The expression for the area in (]72f) looks 
like an action from classical mechanics. If we define the Lagrange density £.: 




then the Euler-Lagrange equation is ( |Goldstein, 1980| ) 



d dC d dC dC 

+ T^TvTTTY = ( 74 ) 



dx d(d x h) dy d(d y h) dh ' 
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FIG. 5 The original minimal surface M and the comparison surface M. 



or, by inserting ([T3|) 







d T h 



+ 



d 



d y h 



0. 



(75) 



dx \^/TTJyW J d v Vv /1 + ( v/i ) 2 , 

This equation has an important interpretation in terms of the surface normal. If we construct 
the two tangent vectors to the surface as we did above, we have 



ui = <9 Sl X(si,s 2 ) = [1,0,^/i] 
u 2 = <9 S2 X(si, s 2 ) = [0, 1, dyh] 



(76) 



and so the surface normal is parallel to Ui x u 2 = [— d x h, —d y h, 1]. Normalizing this vector 
gives us n(si, s 2 ) = — [d x h, d y h, —1] / ^/l + (V/i) 2 and so (|75| ) becomes V-n = 0. In other 
words, a surface which extremizes its area has a divergence-free unit normal. We will address 
this further in the following section. Such a surface is called a minimal surface. Typically, 
(I7SP is written equivalently as ( [Nitschc, 1989| ) the minimal surface equation: 



[1 + {dyhf] dlh ~ 2d X hdyhd X dyh + [l + ^hf] fij/l = 0. 



(77) 



Though we were seeking a surface with a minimum area, the Euler-Lagrange equations 
only guarantee an extremum - we might have a saddle-point or, worse, a maximum! Rigorous 
proofs that the area is minimized in general are difficult. However, we can prove a restricted 
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result easily. If the projection of the boundary DM onto the xy-plane is convex, then it 
is clear that any overhangs on the surface add extra area and that the surface need not 
"stick out" past the boundary of M projected onto the xy-plane. Thus we can use the 
height representation and solve (|77|) subject to a boundary condition h(x,y) = f(x,y) on 
the boundary in the xy-p\a,ne. As depicted in Figure 5, we consider an arbitrary surface M 
with the same boundary as the surface of which we are interested. Since the same range of 
x and y map out both M and M, we can define n(x, y) on M and in the volume between 
the surfaces. Define 

T= / n(x)-dS (78) 
Jm 

where n(x) is the unit normal of the original surface M, and dS is the vector element of 
surface area, pointing along the normal to the surface n(x). If we consider the same integral 
for the original surface M, then since n(x) is parallel to dS, we have n(x) • dS = dS the 
element of surface area! Integrating this over the whole surface just gives us the area A. 
Thus 

A= [ dS= [ n(x)-dS. (79) 
Jm Jm 

Moreover, we have the inequality n(x)-n(x) < 1, since the dot-product of two unit vectors 
is the cosine of some angle and is therefore less than or equal to 1. We arrive at: 

T= f n(x)-dS= / n(x)-S(x)ci£< [ dS = A, (80) 
Jm Jm Jm 

where A is the area of the comparison surface. We are almost done. Note that if we construct 
a closed surface M by gluing together M and M along their common boundary, then the 
integral of n(x) over the whole surface is: 

/ n(x)-dS= / n(x)-dS- [ n(x)-dS = A-T (81) 
Jm Jm Jm 

where the relative minus sign arises because we want the normal to the closed surface to 
always point outward. But by Gauss's law, the first integral in (]81[) can be converted into 
an integral over the volume V enclosed by M, so 

A-T= I n(x)-dS= I V-n(x)dy = 0, (82) 
Jm Jv 
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FIG. 6 Edge (left) and screw (right) dislocations in a smectic- ^4 liquid crystal. 

since M extremizes the area so V-n(x) = 0. Using the inequality (|80| ) we have A = T < A 
and thus the comparison surface has an area at least as large as the original surface M. 
Since M was arbitrary we have shown that our surface M minimizes the area 3 . 



Examples of minimal surfaces abound. The smectic- A liquid crystal phase is a one- 
dimensional crystal, composed of periodic layers of spacing d. Each layer is fluid-like, with 
no internal ordering. Thus, in some sense, each layer behaves as a minimal surface and the 
layer spacing is set by a compression modulus QKamien and Lubensky, 1999|) . The smectic 



phase can have topological defects of two types. As shown in Figure 6, an edge dislocation 
adds or ends a layer in the bulk, while a screw dislocation connects adjacent layers together. 
A screw defect can be represented by the height function: 

h(x,y) = ^tajr 1 (^) (83) 



2tt \x, 

where n is an integer. By choosing n to be an integer, we are guaranteed that when going 
around the origin once we move up exactly n layers. It is straightforward to check that the 



screw dislocation is a minimal surface by direct calculation of (\T7\) . 

Another nice example also arises in smectic liquid crystals. Scherk's first surface (shown 
in Figure 7) QScherk, 1835 ) is a minimal surface which smoothly connects two smectic regions 



which are rotated with respect to each other. It is composed of an infinite set of parallel, 
screw- like dislocations ( |Gido et a/., 1993 ; Kamicn, 2001 ) which must be slightly squashed to 



make the surface satisfy (|77|) . For those who are interested, the height function is defined 



In this context, the vector field n(x) is called a calibration. The interested (or disinterested) reader will 



find this and other results in (Morgan, 2000) 
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FIG. 7 Scherk's first surface. This minimal surface connects together two layered structures with 
different orientations. 

for arbitrary rotation angle a: 



It is interesting to note that as a — ► the height function becomes that of a simple screw 
dislocation fl83|) with nd = 2n. 

B. Mean and Gaussian Curvature: Energetics of Membranes 

When we discussed curves in three dimensions there were not many options for the 
definition of curvature: there was only one coordinate on the curve, and if we wanted 
something that was independent of our parameterization, we need only take derivatives 
with respect to arclength. As we have seen in the last section, however, surfaces can be 
reparameterized and have many more degrees of freedom. As a result, there is more than 
one way to define curvature, in fact there are two. 

We have already hinted at one definition in the discussion of minimal surfaces. For curves, 
we know that a straight line minimizes the distance between two points in space, or in other 
words, k(s) = is the equation of a length-minimizing curve. Generalizing this to surfaces, 
we might define one curvature which is for surfaces that minimize the area spanned by a 
given boundary. Indeed, this is called the mean curvature, H and H = — |V-n(x) where 
n(x) is the surface normal. The other type of curvature is the Gaussian curvature, K, which 
we will define in the following. 




(84) 
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In order to understand this better and to introduce K, we need to step back and consider 
our surface X(si, s 2 ). Since the curvature of a line is k(s) = — T(s) • N'(s), we would like to 
take derivatives of the surface normal along tangent directions of the surface. As in §IIIB, 
we now reintroduce the basis vectors perpendicular to n(si, s 2 ), e 1 (s 1 , s 2 ) and e 2 («i, s 2 ), of 
unit length and mutually orthogonal. We can take directional derivatives 4 of n(si, s 2 ) along 
Gi(s\, s 2 ), (e* • V). Now instead of one number, we have four, which form a matrix: 



ei • [e 2 • V] n e 2 • [e 2 • V] n 

The two dot products in the matrix are somewhat confusing. To be concrete we write this 
matrix using its indices: 



This matrix is also a tensor, known as the Weingarten Map or second fundamental form 5 . 
We can diagonalize this matrix via a similarity transform A = S _1 LS. In this diagonal 
basis the entries of A are precisely what we are after: the first entry in the upper left is the 
derivative of the surface normal along a direction ei. Moreover, since the upper right entry 
vanishes, this derivative has no components along e 2 . Therefore the upper left corner is the 
curvature of the curve in the surface at (s±, s 2 ), tangent to ei. Similarly, the lower right entry 
of A gives us another curvature. Note that the two directions ei and e 2 remain orthogonal 6 . 
They are known as the principal directions on the surface and their associated curvatures 
Ki and k 2 are the principal curvatures (see Figure 8). Equivalently, we can define the two 
principal radii of curvature through R4 = We can easily extract these curvatures from 
the original matrix L. Note that k±k 2 = det A and K\ + k 2 — Tr A. Moreover, since the 

4 It may seem odd to use the gradient in these definitions since n(si, S2) is a function of s\ and S2- The 
derivatives of n(si,S2) along the surface tangent vectors 9 Si X(si,s 2 ) arc just d Sj n(si, s 2 ). By writing 
e»(si,S2) as the linear combination ei = Aijd Sj X., the chain rule implies that (e* • V) = Aijd Sj . 

5 The metric tensor g is sometimes called the first fundamental form. 

6 The orthogonality of the two vectors ei and e 2 follows from the fact the S _1 = S T . In other 
words, L can be diagonalized via an orthogonal transformation since it is symmetric: L12 — £21 = 



ei ■ [ei • V] n e 2 • [ei • V] n 



L = - 



(85) 




(86) 




a 



= n • [V x n] = when n is normal to a surface, i.e. n = V<£/|V</>|. 
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FIG. 8 Curves on a surface (heavy lines) with tangents along the principal directions at the point 
of intersection. 



trace is cyclic and det AB = det A det B, we see that: 

Tr L = K 1 (s 1 ,s 2 ) +k 2 (si,s 2 ) 

detL = «i(si, s 2 ) ac 2 (si, s 2 ). (87) 

The product of the curvatures is known as the Gaussian curvature, K = K\K 2 , while the 
average of the curvatures is the mean curvature, H = | [k\ + k 2 ] . The Gaussian and mean 
curvature contain all the information to describe the bending of our surface. 

The reader may be concerned that our definitions of curvature remain dependent on our 
choice of basis vectors. This is not the case. Indeed, using the expression (|86D we have 

2 



H = ~> efef 



dn 



a 



2 ^ 1 1 dxP 
i=i 

K = (e^elei - e^ 2 ele{) §^ (88) 

These expressions still look dependent on ej. However, since n is a unit vector, n a rfid^n a = 
and so we have 



2 



1 dn a 1 ^ n , n dn a 1 



where we have used the fact that {ei, e 2 , n} form an orthonormal triad. Thus, as we hinted 
at the beginning of this section, the mean curvature is proportional to the divergence of the 
surface normal. 
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We have to do a little more work on the expression for the Gaussian curvature. The form 
in ( |38D reminds us of the expressions in the derivation of the Mermin-Ho relation. Using the 
orthonormality of the basis vectors, 

K = e a x el [e{e s 2 - efef) d^dfli 1 = e^e1e^ p n p dpn a d 5 n^ '. (90) 

Finally, as in the discussion of (fiO"D, we have 

K = ^e aPl n~<e ijk n k d i n a d j n f3 = n • [V x ft] . (91) 

Moreover, the Gaussian curvature has a simple interpretation. Consider the normal spherical 
map defined by analogy with the tangent spherical map. For each point on the surface we 
identify a point on the unit sphere which corresponds to the surface normal at that point. 
This map is also known as the Gauss map. Since the surface element on our surface is just 
dS/j, = n^dS, we see from the discussion at the end of §IIIC, that the Gaussian curvature 
is just the area swept out by the surface normal on the Gauss map as we move along the 
surface. In other words, if we have a small region M of our surface, then 

/ KdS= I -€ a p 1 ri 1 € ijk d i n cl dj'nP dS k . (92) 

J M JM * 

Thus the Gaussian curvature is the ratio between the infinitesimal area swept out on the 
Gauss map and the infinitesimal area of the original surface to which it corresponds. 

We have seen that the two principal curvatures depend only on the surface normal and not 
on our choice of coordinates or basis vectors. This is useful when we consider the energetics of 
fluid membranes. They are called "fluid" because they have no internal structure. Therefore 
it would be unphysical to build an energy out of anything but the two invariants H and K. 
Note that if we have an open surface there is no distinction between inside and outside, so 
the layer normal n(s\,S2) is defined only up to a sign. Though the Gaussian curvature is 
independent of this sign, the mean curvature is not. To get around this, we insist that the 
free energy be even in powers of H flCanham, 1970| ; [Helfrich, 1973|) : 



Fqh = J dS {2kH 2 + kK} , (93) 

where k and K are (confusingly) the standard symbols for the bending moduli. This is known 
as the Canham- Helfrich free energy for fluid membranes. Note that they both have units of 
energy, since the dimensions of K and if 2 cancel the dimensions of the surface. Also note that 
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the integration is done with respect to the actual surface area, so that dS = y/det gds\ds2- 
We will see in the following sections that the integral of the Gaussian curvature is a constant 
and so the term proportional to ~k is usually neglected. We will also discover that not only 
is the Gaussian curvature independent of our choice of basis vectors, but that it also can 
be measured with no knowledge of the surface normal n(sj, S2). Finally, if the membrane is 
tethered then it has internal elastic degrees of freedom. These degrees of freedom couple to 
the geometry and produce a variety of singular structures ( piDonna et a/., 2002| ). 

V. GLOBAL THEORY OF SURFACES 

A. The Gauss-Bonnet Theorem: Foams on Curved Surfaces 

We have seen that there are two sorts of curvature that we can consider on a surface, the 
mean and Gaussian curvatures. We have studied both of these by considering deviations of 
the normal vector to the surface. However, a large part of differential geometry focuses on 
intrinsic properties, those quantities which can be measured without reference to the space 
in which the manifold is embedded. It turns out that the Gaussian curvature is intrinsic: 
by measuring the lengths and diameters of small circles entirely in the surface, one can 
determine K. 

We have all the technology necessary to demonstrate this remarkable property. To show 
this, we consider a patch on our surface M, with boundary dM. As in §IVA, we can take 
the unit normal at each point and map it to the unit sphere. As shown in Figure 9, when 
we traverse the patch M on the surface, we traverse a patch M on the unit sphere. As we 
discussed, the area of M is the integral of the Gaussian curvature of the patch M: 

ffdS = [[ KdS= [[ \e^ p e aPl n a ^)d v n\^)d p n\^)dS^ (94) 
JJm JJm JJm * 

where we have used (|9T|). Now consider the boundary curve dM. This curve has a tangent 

vector T(s) in space. We know from our discussion of curves that T'(s) = «(s)N(s). 

However, while T(s) is also tangent to the sphere, N(s) is, in general, not. The normal to the 

curve, N(s), can have a component along the normal to the surface n(s) and perpendicular 

to it. If we lived in the surface and could only make measurements in the surface, like some 

kind of bug, then we could only measure the component of the curve's normal in the surface, 
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FIG. 9 Geometry for the Gauss-Bonnet Theorem. Note that now we are using the Gauss map 
which maps the surface normal n(x) unto the unit sphere. 

the surface normal "y(s): 

7(S ) = T-M [T- M .n M ln( S ) (g5) 
KW^/l-INW-nW] 2 

This is a unit vector which lies in the plane tangent to the surface and which is perpendicular 
to T(s), that is, the tangent to the boundary dM of the patch M. We would like to define 
a curvature which measures the extra curvature in the curve, not arising from the curvature 
of the surface in which it is embedded. By analogy to the Frenet-Serret formulae, we define 
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the geodesic curvature 7 n g to be 

k 9 {s) = T(s) • 7 ( a ) = K (s)y/l-\N(8)-n(s)] 2 . (96) 

Since {T(s), 7(5), n(s)} form an orthonormal triad, we have 

k 9 (s) = T'(s) • [n(a) x T(s)} = n(a) ■ [T( S ) x T'( S )] , (97) 

a form we have seen before in the context of the Mermin-Ho relation! Since the unit tangent 
vector lies in the plane perpendicular to n(x), we may write it in terms of our arbitrary 
basis ei(x) and e 2 (x), introduced in §IIIB, T(s) = cos^ex + sin6 , e 2 . Using (0) to calculate 
Kg, we find that 

k 9 (s) = d s 6(s) - e?(s)d a e%(s) (98) 

where we have used the fact that for a unit vector u(s), d s [u(s)] 2 = 2u(s) • d s u(s) = 0. In 
analogy with (El), we can integrate the geodesic curvature around the boundary of M to 



find: 



/ /t 9 (s)rfs = <f [V0(x) - n(x)] • dR. (99) 

As with our discussion of the relation between link, twist and writhe, we may transform the 
line integral of f2(x) into a surface integral over M so that 

I d s 9{s)ds=l K g (s)ds+ [[ KdS (100) 

JdM JdM JJm 



where we have used ([IT]) to rewrite the final integrand as the Gaussian curvature 
Finally, if the boundary curve does not intersect itself, then the tangent vector rotates around 
n(x) exactly once, so 9(s) changes by 2n around the curve. We have thus established the 
Gauss-Bonnet Theorem: 



Kds + <b Kgds = 2tt. (101) 

M JdM 

By integrating the geodesic curvature around a closed loop, we can calculate the integrated 
Gaussian curvature that we surround. As the loop shrinks ever smaller around a point Xq, 



7 We may also define the normal curvature, K n (s) which is the curvature of our embedded curve that is 
imposed by the surface. We define K n (s) — T'(s) ■ n(s) and thus the total curvature satisfies n(s) — 

^Kl(s) +K 2 g (s). 
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we can calculate the Gaussian curvature at that point through division by the surface area 
enclosed. This is remarkable since the geodesic curvature can be measured without the use of 
the layer normall In other words, since the geodesic curvature is intrinsic, so is the Gaussian 
curvature! We can take this result a step further by considering regions that have discrete 
angles in their boundaries (i.e. polygons). When there is a sharp bend in the boundary 
curve dM, we cannot calculate the geodesic curvature k 9 (s). We can, however, integrate 
the geodesic curvature along the smooth parts of the boundary. If there are j sharp bends 
with angles A9j, respectively, then when we integrate around the boundary as in (|99"D we 
will have a deficit of A9j, or in other words the smooth part of 6(s) does not have to 
change by 2n for the curve to come around. We thus have 

/ f KdS + <f K g ds + V A9j = 2tt (102) 

JjM JdM ■ 

where we understand that the integral around the boundary should be broken into smooth 
segments of the boundary. The jump angles account for the discontinuities. Note that we 
can have lines for which K g (s) = 0. These are the "straight" lines on the surface and are 
called geodesies. If we were to build a polygon with n-sides, all of which are geodesies, then 



SM) reads: 



2tt-V A8, = // KdS. (103) 



j - JM 



Since the sum of the external angles A9j of a polygon in flat space is 2n, we see that the 
Gaussian curvature is a measure of the excess (or deficit) angle in a polygon. If the Gaussian 
curvature is positive the sum of the external angles is smaller than 2tt so the sum of the 
internal angles Ylji 71 ~ ^j) = ( n ~ 2)7T + jjKdS is larger than we might expect. For 
example, imagine the triangle on the sphere connecting the North Pole (90°N), Pontianak, 
Indonesia (roughly 0°N, 109°20'W) and Loolo, The Congo (roughly 0°S, 19°20'W) along 
great circles. At each vertex of this triangle the arcs of great circles meet at 90° and so the 
sum of the interior angles is 37r/2. Since n = 3 for a triangle, we discover that jj KdS = it/ 2 
for this triangle that covers ~ of the globe. However, we know that for a sphere of radius R, 
the Gaussian curvature is just 1/R 2 and so jj KdS = |^r- = vr/2, as we expect! 

There are many beautiful uses and examples of the interplay between Gaussian curva- 
ture and geodesic curvature in physical systems. For instance, Avron and Levine (Avron 



and Levine, 1992]) have considered dry foams on curved, two-dimensional surfaces. In this 
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context, "dry" refers to the fact that there is no fluid between adjacent bubbles so that 
the walls between them can be treated as lines and the vertices may be treated as points. 
There are two key ingredients to the physics of dry foams: surface tension and pressure. In 
two dimensions the surface tension amounts to a line tension o along the interfaces between 
neighboring bubbles. If the soap film is uniform then the line tension is constant as well. 
Whenever the lines meet, mechanical equilibrium must be maintained. The pressure in the 
bubbles exerts a force on their boundaries. If two cells are separated by a boundary line then 
the pressure difference AP must be balanced by the boundary. As with the derivation of the 
wave equation, the force exerted by the boundary curve T is aT'(s). However, the pressure 
exerts forces only within the surface and so we are only interested in the component of the 
force in the plane of the surface. According to the above discussion, that force is aK g (s)~y(s), 
and so the magnitude of the force is <jk 9 (s) = AP. In the case of planar foams, K g (s) = k(s) 
and this is known as the Young-Laplace law. To model the diffusion of gas from one cell to 
the other, we assume a simple dynamics where 

^ = -C$>P/, (104) 

3 

where N is the number of gas molecules in the cell of interest, APj is the pressure difference 
between it and its j th neighbor, ij is the length of the boundary separating the cell from 
its j th neighbor, and C > is a diffusion constant. This dynamics captures the simple idea 
that if a bubble is higher pressure than its neighbors so that APj > 0, then it loses gas, 
while if it has lower pressure it gains gas. Using the generalization of the Young-Laplace 
law and the Gauss-Bonnet theorem we have: 



^L = -Ca(f K g (s)ds = Ca\ [[ KdS + Vl 

dt JdM yJjM ~* 



7t - aj ) - 2tt } (105) 



where aj are the internal angles. 

For flat surfaces with K = it is known that a hexagonal honeycomb network of bound- 
aries minimizes the length of the cell walls ( [Hales, 200 1 ; [Morgan, 2000| ). If we imagine 
starting with a stable configuration on a flat membrane and distorting the membrane, we 
should take each internal angle to be 27r/3. Doing so Avron and Levine ( |Avron and Levine" 



T992|) found 

c 

~dt 



= C<r{~{n-6) + I/ KdS}. (106) 

M 
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If we seek a stationary configuration so that N is time independent, we see that for flat 
membranes K — and only bubbles with n = 6 sides are stationary: hexagons. If the 
surface has positive curvature, K > 0, and there is an instability. To make (dN/dt) = 0, we 
see that n < 6. However, even if an area and an n were found to make N time independent, 
we can see that there is an instability: if the area grows then (dN/dt) becomes positive so 
the bubble grows more. Likewise, if the bubble shrinks then the integral over K gets smaller 
and so gas flows out of the bubble and it shrinks some more. Thus on a positively curved 
surface the only stable situation is one for which one bubble overtakes the whole system. On 
the other hand, when the surface has negative Gaussian curvature, it follows from a similar 
argument that when the area of the bubble increases (dN/ dt) < and when the area shrinks 
(dN/dt) > 0. In the critical case of a flat membrane the stationary solution is a hexagon 
and it is neither stable nor unstable towards growth ([Avron and Levine, 1992| ). 



B. The Euler Characteristic and the Genus: Defects on Surfaces 



We can take (102) one step further by considering a closed surface. In this case we 



can integrate the Gaussian curvature over the whole manifold. If we triangulate the entire 
surface to form a net, then we can use the Gauss-Bonnet theorem to establish a relation 
between the topology of the network and the total Gaussian curvature. At each vertex there 
is a total angle of 2n which is divided into the internal angles of the triangles meeting at that 
point. Each face contributes 7r to the total angle at the vertices in addition to the excess 
angle from that face, JJ KdS. Adding all the triangles together we have: 



2ttV = ixF KdS (107) 

where V is the number of vertices and F is the number of faces. Each face contributes three 
edges E, but each edge is shared by two triangles, so 3F = 2E. We then have 

V - E + F = —(ft KdS. (108) 
2?U7 m 

This is a remarkable result. As we show in the bottom two graphs of Figure 10 , if we remove 
an edge from our network (E — ► E — 1), two faces join into one (F — > F — 1), and we lose two 
vertices (V — > V — 2), so V — E + F does not change! Likewise if we add an edge V — E + F 
is unchanged. Thus our result does not require the use of triangles. This invariant is known 
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as the Euler characteristic \ = V — E + F. Note that if we have any network that can be 
deformed into a sphere without cutting the edges or changing the vertices, then \ can be 
calculated using ( |108| ) for a perfect sphere, where the surface area is 4nR 2 and the Gaussian 
curvature is 1/R 2 . Thus we find that 

V-E + F = X = 2 (109) 

for a network with the topology of a sphere. Suppose we take two faces on the sphere, 
deform them to be the same triangle and place them together, face-to-face. In the process 
we lose three edges (E — > E — 3), two faces (F — > F — 2) and three vertices (V — > V — 3) 
so x = 0- What have we done? We have made a doughnut (or a torus) with one handle 
and have reduced the Euler characteristic by two. Clearly, any time we add a handle x is 
reduced by two. If we define the genus g to be the number of handles of the surface, then 

1-^KdS = V-E + F = X = 2(l-g). (110) 

Though the Gaussian curvature was a geometric property, when integrated over the entire 
surface it becomes a topological invariant, independent of the local geometry. 

The Euler characteristic can be used to understand the topology of defects on closed 
surfaces. Suppose that we have a unit vector field v(x) living in the local tangent plane 
to a closed surface of genus g. Since v(x) lies in the tangent plane, we have v(x) = 
cos#(x)ei(x) + sin#(x)e2(x), where we have reintroduced our vectors ei(x) and e2(x) which 
are everywhere perpendicular to the unit normal n(x) of the surface. We will try to cover the 
surface with a vector field that is single valued. In other words, if we integrate derivatives 
of v(x) around a closed curve, we should get back the same vector. But as we discussed 
in §IIIB, to ensure this property we should focus on the covariant derivative D0(x). We 
will have a single-valued vector field when the curl of this derivative vanishes. In this case 
equations ( pH]) and (|9"ID give: 

[VxV0(x)]-rfS = tf KdS = 2ttx (HI) 

M JJ M 

Usually V x V# = 0, but this result tells us that it is not true on an arbitrary genus surface 
M. 
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FIG. 10 Removing or adding an edge from a graph maintains the value of the Euler characteristic 
V -E + F. 




T 

FIG. 11 A vector field on a triangular patch. Note that the +1 defects at the vertices and at the 
center force —1 defects on the edges. 
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Without presenting a detailed discussion on topological defects flChaikin and Lubensky 



1995| ), we can see that V x V0(x) does not vanish in the presence of a defect configuration 



and that on a two dimensional surface: 

VxV«(x)=^m/(x-Xi), (112) 

i 

where rrii is the charge of the defect and Xj is the corresponding position. To see this, we 
consider 

0(x) = mtan" 1 (|) . (113) 

It is straightforward to calculate V#(x) and then to integrate the gradient around a closed 
curve 7 that contains the origin. We have 



2nm = <p V0(x) • dx 

V x V0(x) • dS, (114) 



7 



'M 

where 7 is the boundary of M. Since we may arbitrarily shrink the path 7 around the origin 
and the integral remains constant, it must be that V x V0 = m5 2 (x). Generalizing this to 
an arbitrary collection of defects leads us to (|112|) . Using (|112|) in ( |111| ) gives us: 

£mj = 2(1 -(/) = *. (115) 

i 

This result is known as the Poincare-Brouwer theorem. Since a vector field must be single 
valued on the surface, rrii must be an integer so that cos 9 and sin 9 are well-defined. This 
result tells us, for instance, that a vector field on the surface of a sphere must have two +1 
defects or one +2 defect. On a torus, however, no defects are necessary. 

There is another, more geometric way to establish ( |115|) . We first argue that two vector 
fields on the same surface must have the same total topological charge. Consider two vector 
fields u(x) and v(x). We can triangulate the surface so that each triangle contains at most 
one defect in u(x) and one defect of v(x) and that every defect is in some triangle. In each 
triangle we may now integrate V#„(x) and V0„(x) around the triangle boundary. The total 
topological charge \u of u(x) is 



E E / V0„(x)-dx (116) 

: 1 — „,J~~~ ^ 



Xu 

triangles edges 
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FIG. 12 A vesicle with some sort of vector order parameter. The Poincare-Brouwer theorem 
assures that there are two defects in the vector field, drawn as the dark lines on the vesicle. The 
vesicle distorts so that the Gaussian curvature is larger in the vicinity of the defects. 

and similarly for v(x). Note that to establish the angles 6 u (x) and #„(x), we must choose 
basis vectors e^x). Since we know that V# u (x) depends on our choice of these vectors, the 
integrals of V# u (x) around each triangle depend on the underlying geometry through the 
Mermin-Ho relation ([Hp. More importantly, each vector field on M may require a different 
choice of basis vectors e$ (x) , since a defect in any particular triangle forces a constraint on 
ej(x). Thus it is not, in general, possible to calculate Xu and Xv using the same set of ej(x). 
However, the difference between Xu and Xv is just 



The difference vanishes since the angle between u(x) and v(x) is independent of our choice 
of basis vectors. Therefore when we sum over all triangles we get each edge twice but in 
opposite directions. Since we do not need the basis vectors to define the integrals, it is clear 
that the integrals in opposite directions cancel and we get 0. Thus any two vector fields 
have the same topological charge. 

We will now count the defects of a vector field that we construct on triangulated surface, 
via the following rules: (1) put a defect with +1 charge at each vertex, (2) put a defect 
with +1 charge at each center, and (3) put a defect with —1 charge on each edge. We can 
see that this will produce a consistent vector field everywhere else, as shown in Figure 11. 




triangles edges 



(117) 
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Adding together the charges we have 



Xu = V-E + F = X , (118) 

and so the topological charge is the Euler characteristic for any vector field on the surface. 

A striking example of the effect of these defects can be seen when there is a coupling 
between the bending of the surface and the configuration of the vector field (|MacKintosh 



and Lubensky, 1991| ; park et a/., 1992|) . For instance, if we have a vector field living on a 
closed surface then it must have two +1 defects. But because the defects are discrete, when 
we consider (|1 1 1|) we find that the Gaussian curvature K must also be confined to the core 
of the defect since we may also integrate fl4~l|) over any submanifold of M. But then all the 
Gaussian curvature sits at the defects. If instead of viewing (41) as a constraint, we write 
it as an energy: 

F = \j dSA[V9(x) - fi(x)] 2 + F CH , (119) 



where Fqk is the free energy for the surface fl93|), and A is the spin-stiffness. On a flat 
surface we can always choose Q = and so the first term in ( |119| ) just describes the usual 
Goldstone mode of broken rotational invariance. On a curved surface the inclusion of Q is 
necessary in order to make the free energy uniquely defined. The energy ( |119j ) now represents 
a competition between surface bending and the requirement of two +1 topological defects 
which serve as a source for Gaussian curvature. Both of these energies vanish for the plane, 
but if we restrict our topology to a sphere, for instance, then neither term can be made to 
vanish everywhere. Balancing these two effects, a vesicle would form an oblate shape, as 
shown in Figure 12. 



VI. THREE DIMENSIONS AND BEYOND 



We have covered the basic elements of the geometry of curves and surfaces. Things become 
significantly more abstract in three dimensions because there is no analog of the normal 
vector, binormal vector or surface normal: the three dimensional system is not embedded in 
a higher dimensional space. Though it might be hard to visualize how something can have 
curvature without an embedding space, we already know that it is possible. Recall that 
the Gaussian curvature could be measured without any knowledge of the surface normal. 
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Quantities with that property are called intrinsic. The general theory of relativity relies on 
intrinsic quantities to describe the curvature of four-dimensional space-time QMisner et at., 



|1973| ) - there is no reference to a larger space in which our universe lives. Though we will 
not go into any mathematical detail, it is worthwhile to describe two examples from soft 
materials. 

The first is the blue phases of chiral liquid crystals ( [Meiboom et at., 1981fc |Wright and 



Mermin, 1989Q . These are phases in which there is a three-dimensional, periodic modulation 
of the nematic director n(x) with a lengthscale comparable to that of visible light. These 
usually exist only over a narrow temperature range and are stabilized by an often neglected 
term in the nematic free energy. The Frank free energy for a nematic liquid crystal is 

F = ^ J d 3 x (V • n) 2 + K 2 [n ■ (V x n)] 2 + K 3 [n x (V x n)] 2 

+2K 2i V • [(n • V) n - (V • n) n] } (120) 

The standard elastic constants Ki, K 2 and K 3 are a measure of the energy cost for splay, 
twist and bend modes, respectively ( |Chaikin and Lubensky, 1995| ). The last term with 



elastic constant K 24 is known as the saddle-splay term. Though it is a total derivative, 
when there are defects present it can contribute to the energy. The blue phase is riddled 
with precisely those defects that contribute to saddle-splay. The remarkable thing about 
this term is that it is precisely the Gaussian curvature as in (OT): 

n(x) • V x n(x) = — V ■ [(n • V) n - (V ■ n) n] (121) 

What does this mean? If there are surfaces to which the nematic director is normal, then 
the saddle-splay is the Gaussian curvature of those surfaces. However, the saddle-splay is 
more general. Even if the director is not a field of layer normals, the saddle-splay is a 
measure of curvature. Indeed, the blue phases can be understood as unfrustrated systems 
in curved three-dimensional space ( [Sethna et al, 1983 ). Projecting this texture into flat 



space leads to the topological defects in these phases. A compelling and alternative way 
of viewing the blue phases is to view them as decorations of space filling minimal surfaces 
( Pansu and Dubois- Violette, 19891) . It is the connection between topological defects and 
curvature that makes it possible to identify and locate defects in lattice simulations of liquid 
crystals flPriezjev and Pelcovits, 200 1| ). 
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Another three-dimensional system that can be understood in terms of curvature is the 
melting and freezing of hard-sphere fluids. In two dimensions, as hard-spheres condense from 
the fluid phase to the crystalline phase, they form close-packed triangles. Eventually these 
can assemble into a triangular lattice. The situation is not so happy in three-dimensions. 
There, four spheres can close-pack into a tetrahedron, but tetrahedra cannot assemble to fill 
space. As a result, all lattices have a packing density lower than the best local packing. This 
difficulty can be viewed as a geometric frustration along the same lines as the blue phase. 
In positively curved space, however, this frustration can also be eliminated ( [Nelson, 1983 ; 
Nelson and Widom, 1983|) , and tetrahedral close-packing can fill space. The BCC phase is 
the simplest example of this structure in flat space. 



Outlook 



Though the mathematics described here has a certain elegance and beauty, I hope that 
I have conveyed the utility of differential geometry in a variety of physical problems. While 
a geometric description of a system often leads to an intuitive perspective, there are many 
other arenas in which a geometric formulation of the problem is not only useful but essential. 
Put together with statistical mechanics, differential geometry has been and will continue to 
be a powerful tool in the study of soft materials. 



Acknowledgments 

I thank I. Bluestein, B.A. DiDonna, G.M. Grason, M.I. Kamien, W.Y. Kung, T.C Luben- 
sky, T.R. Powers, P. Ziherl and especially G. Jungman and P. Nelson for careful and critical 
reading of this review. This work was supported by NSF Grants DMRO 1-29804 and INT99- 
10017, and a gift from L.J. Bernstein. 



43 



References 



Aldinger, J., I. Klapper, and M. Tabor, 1995, J. Knot Theory Ramifications 4, 343. 
Avron, J. E., and D. Levine, 1992, Phys. Rev. Lett. 69, 208. 
Bensimon, D., D. Dohmi, and M. Mezard, 1998, Europhys. Lett. 42, 97. 
Canham, P., 1970, J. Theor. Biol. 26, 61. 

Chaikin, P. M., and T. C. Lubensky, 1995, Principles of Condensed Matter Physics (Cambridge 

University Press, Cambridge). 
David, F., 1989, in Statistical Mechanics of Membranes and Surfaces: Jerusalem Winter School 

for Theoretical Physics, edited by D. R. Nelson, T. Piran, and S. Weinberg (World Scientific, 

Singapore) . 

DiDonna, B. A., T. A. Witten, S. C. Venkataramani, and E. M. Kramer, 2002, Phys. Rev. E 65, 
016603. 

Doi, M., and S. F. Edwards, 1986, The Theory of Polymer Dynamics (Oxford University Press, 
Oxford) . 

Fain, B., and J. Rudnick, 1999, Phys. Rev. E 60, 7239. 

Fain, B., J. Rudnick, and S. Ostlund, 1996, Phys. Rev. E 55, 7364. 

Fuller, F. B., 1971, Proc. Natl. Acad. Sci. USA 68, 815. 

Fuller, F. B., 1978, Proc. Natl. Acad. Sci. USA 75, 3557. 

de Gennes, P.-C, 1970, Scaling Concepts in Polymer Physics (Cornell University Press, Ithaca, 
NY). 

Gido, S. P., J. Gunther, E. L. Thomas, and D. Hoffman, 1993, Macromolecules 26, 4506. 
Goldstein, H., 1980, Classical Mechanics, 2 nd ed. (Addison- Wesley, Reading). 
Goldstein, R. E., T. R. Powers, and C. H. Wiggins, 1998, Phys. Rev. Lett. 80, 5232. 
Hales, T. C, 2001, Discrete Comput. Geom. 25, 1. 

Harris, A. B., R. D. Kamien, and T. C. Lubensky, 1999, Rev. Mod. Phys. 71, 1745. 

Helfrich, W., 1973, Z. Naturforsch 28C, 693. 

Julicher, F., 1994, Phys. Rev. E 49, 2429. 

Kamien, R. D., 1998, Eur. Phys. J. B 1, 1. 

Kamien, R. D., 2001, Appl. Math. Lett. 14, 797. 

Kamien, R. D., and T. C. Lubensky, 1999, Phys. Rev. Lett. 82, 2892. 



44 



Kamien, R. D., T. C. Lubensky, P. Nelson, and C. S. O'Hern, 1997, Europhys. Lett. 38, 237. 
Kleinert, H., 1990, Path Integrals in Quantum Mechanics Statistics and Polymer Physics (World 

Scientific, Singapore), chap. 16. 
MacKintosh, F. C, and T. C. Lubensky, 1991, Phys. Rev. Lett. 67, 1169. 
Marko, J. F., 1997, Europhys. Lett. 38, 183. 
Marko, J. F., and E. D. Siggia, 1994, Macromolecules 24, 981. 
Marko, J. F., and E. D. Siggia, 1995, Phys. Rev. E 52, 2912. 

Meiboom, S., J. P. Sethna, P. W. Anderson, and W. F. Brinkman, 1981, Phys. Rev. Lett. 46, 1216. 
Mermin, N. D., and T. L. Ho, 1976, Phys. Rev. Lett. 36, 594. 

Millman, R. S., and G. D. Parker, 1977, Elements of Differential Geometry (Prentice-Hall, Engle- 
wood Cliffs, NJ). 

Misner, C. W., K. S. Thorne, and J. A. Wheeler, 1973, Gravitation (Freeman, San Francisco). 

Morgan, F., 2000, Geometric Measure Theory, 3rd ed. (Academic Press, San Diego). 

Moroz, J. D., and R. D. Kamien, 1997, Nucl. Phys. B 506, 695. 

Moroz, J. D., and P. Nelson, 1997, Proc. Nat. Acad. Sci. USA 94, 14418. 

Moroz, J. D., and P. Nelson, 1998, Macromolecules 31, 6333. 

Nelson, D. R., 1983, Phys. Rev. B 28, 5515. 

Nelson, D. R., and M. Widom, 1983, Nucl. Phys. B 240, 113. 

Nelson, P., 1998, Phys. Rev. Lett. 80, 5810. 

Nitsche, J. C. C, 1989, Lectures on Minimal Surfaces, (Translated by J.M. Feinberg) (Cambridge 

University Press, Cambridge). 
Pansu, B., and E. Dubois- Violette, 1989, Europhys. Lett. 10, 43. 
Park, J., T. C. Lubensky, and F. C. MacKintosh, 1992, Europhys. Lett. 20, 279. 
Priezjev, N. V., and R. A. Pelcovits, 2001, Phys. Rev. E 64, 031710. 



Scherk, H. F., 1835, J. r. angew. Math. 13, 185, also see flNitsche, 1989[) . 

Sethna, J. P., D. C. Wright, and N. D. Mermin, 1983, Phys. Rev. Lett. 51, 467. 

Vologodskii, A. V., S. D. Levene, K. V. Klenin, M. Frank-Kamenetskii, and N. R. Cozzarelli, 1992, 

J. Mol. Biol. 227, 1224. 
Watson, J. D., and F. H. C. Crick, 1953, Nature 171, 737. 

Wolgemuth, C. W., T. R. Powers, and R. E. Goldstein, 2000, Phys. Rev. Lett. 84, 1623. 
Wright, D. C, and N. D. Mermin, 1989, Rev. Mod. Phys. 61, 385. 

45 



